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PREFACE 


An  understanding  of  the  radiatSfen  signatures  of  reentry  vehicles 
is  of  fundamental  importance  for  both  offensive  and  defensive  missile 
systems.  This  Memorandum  represents  an  extension  of  our  knowledge 
of  reentry  radiation  to  a  new  area,  that  of  nonequilibritmi  wake 
radiation.  This  work  is  part  of  a  continuing  theoretical  study  of 
reentry  physics  supported  by  the  Advanced  Research  Projects  Agency. 


SUMMARY 


The  wakes  of  reentry  vehicles  at  altitudes  of  100,000  ft  or 
above  will  not  be  in  chemical  equilibrium.  Under  these  conditions 
the  wake  radiation  cannot  be  obtained  from  the  ordinary  equilibrium 
calculations,  but  a  detailed  study  must  be  made  of  the  chemical 
kinetics  of  the  nonequilibrium  wake.  This  Memorandum  presents  the 
results  of  such  a  study  for  the  case  of  a  blunt,  nonablating  vehicle. 

In  this  nonequilibrium  regime  the  radiation  (UV,  visible,  and 
near  IR)  will  be  produced  primarily  by  various  chemiluminescent 
reactions.  For  the  case  of  a  pure  air  wake  the  most  important  of 
these  will  be  the  Lewis-Rayleigh  nitrogen  afterglow,  the  blue  nitric 
oxide  afterglow,  and  the  NO^  "continuum."  The  amount  of  radiation 
produced  by  these  reactions  will  depend  on  the  species  concentrations 
in  the  wake,  which  will  in  turn  be  affected  by  the  chemical  reactions 
and  the  amount  of  turbulent  mixing  occurring  in  the  wake. 

Knowledge  of  the  important  chemical  reactions  for  a  pure  air 
wake  is  fairly  complete.  However,  large  uncertainties  exist  with 
respect  to  hypersonic  turbulence.  For  that  reason,  two  recently 
proposed  models  of  the  turbulent  mixing  process  are  used  in  this 
Memorandum;  the  inviscid  random  convection  model,  and  the  homogeneous 
mixing  model.  These  two  models  give  radiant  intensities  for  the 
various  chemiluminescent  reactions  that  differ  as  much  as  several  orders 
of  magnitude.  This  large  difference  is  another  illustration  of  the 
importance  of  increasing  our  knowledge  of  hypersonic  turbulence. 

The  radiation  estimates  presented  here  show  that  the  nonequilibrium 
wake  radiation  in  the  UV,  visible,  and  near  IR  should  dominate  the  gas 
cap  and  surface  radiation  at  altitudes  above  about  100,000  ft. 
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I.  INTROD0CTION 


During  the  past  several  years  there  has  been  a  great  deal  of 
interest  in  the  radiation  emitted  during  the  reentry  of  hypersonic 
vehicles.  There  has  been  extensive  work  on  the  various  sources  of 
radiation,  with  particular  emphasis  on  the  radiation  from  the  shock 
layer  at  the  front  of  the  reentry  body,  and  from  the  body  surface 
Itself.  More;  recently,  it  has  been  realized  that  for  most  reentry 
bodies  the  dominant  source  of  radiation,  over  large  portions  of  the 
reentry  path,  will  be  the  wake.  This  realization  has  stimulated 
extensive  work  on  wake  radiation.  To  date,  however,  the  bulk  of  this 
work  has  been  conceded  with  the  radiation  from  a  wake  in  chemical 
equilibrium. 

Recent  studies  of  the  chemical  kinetics  of  air  and  of  chemically 
reacting  hypersonic  flow  fields  have  shown  that  for  altitudes  above 
about  50,000  ft  the  flow  field  behind  a  reentry  vehicle  will  not  be 
in  chemical  equilibrium; as  the  flow  expands  and  cools  behind  the 
body,  the  atomic  recombination  processes  will  lag  behind  the  flow, 
resulting  in  a  nonequilibrium  chemical  composition.  This  nonequilibrium 
chemistry  will  affect  the  temperature  and  density  properties  of  the 
flow  field.  It  will  also  affect  the  nature  of  the  radiation.  For 
slender  bodies,  where  the  initial  degree  of  dissociation  in  the  flow 
is  small,  the  difference  between  the  equilibrium  and  nonequilibriim 
chemistry  will  not  be  very  large,  and  the  resulting  effects  on  the 
flow  field  and  on  the  radiation  will  be  slight.  For  blunt  bodies 
however,  where  the  air  behind  the  bow  shock  is  almost  completely 
dissociated,  the  difference  between  equilibrium  and  nonequilibrium 
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chemistry  will  be  very  large;  it  will  greatly  affect  the  temperature 
of  the  flow  field,  ar^,  it  will  completely  change  the  nature  of  the 
radiation. 

The  purpose  of  this  Memorandum  is  to  begin  the  study  of  radia¬ 
tion  from  nonequilibrium  wakes.  For  a  complete  study  of  this  problem, 
one  should  include  the  effects  of  ablation  products  in  the  wake,  since 
they  can  be  a  very  important  source  of  radiation.  In  this  Memorandum 
ablation  products  will  be  canitted  and  only  pure  air  wakes  will  be 
considered.  The  results  of  this  study,  besides  being  directly  applicable 
to  nonablating  vehicles,  are  also  a  necessary  preliminary  to  a  study  of 
ablation  radiation,  since  before  one  can  hope  to  understand  the  chemistry 
of  the  ablation  product--air  mixture,  one  must  understand  the  pure  air 
problem. 

The  primary  aim  of  this  Memorandum  will  be  to  describe  the 
important  aerodynamic  and  chemical  processes  affecting  radiation  in  the 
nonequilibrium  regime,  to  discuss  the  connection  between  the  flow 
field  properties  and  the  radiation,  and  to  describe  the  characteristics 
of  the  radiation  that  is  emitted.  Numerical  estimates  of  the  intensity 
of  the  nonequilibrium  radiation  will  be  made  for  various  values  of  the 
flight  parameters.  The  only  case  to  be  considered  will  be  that  of  blunt 
bodies,  where  the  nonequilibrium  effects  will  be  the  most  pronounced,  and 
only  radiation  in  the  UV,  visible,  and  near  IR,  will  be  considered,  since 
this  is  the  spectral  region  containing  the  bulk  of  the  radiant  energy. 

Before  beginning  the  detailed  discussion,  it  is  perhaps  useful 
to  summarize  the  main  differences  between  the  equilibrium  and  non- 
equilibritim  situations.  In  the  case  of  an  equilibrium  wake  the 
expansion  of  the  flow  behind  the  body  is  accompanied  by  atomic 
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recombination  with  its  attendant  release  of  energy.  Because  of  this 
energy  release,  the  far  wake  (i.e.,  the  wake  downstream  of  the  P 
point)  will  be  fairly  hot,  with  temperatures  initially  of  the  order 
of  3000°K  to  4000°K.  Hie  far  wake  will  also  be  primarily  molecular;  the 
only  atomic  species  present  in  significant  quantities  will  be  oxygen. 
Because  the  wake  is  assumed  to  be  in  equilibrium,  the  excited  states 
of  the  air  molecules  will  have  an  equilibrium  population  distribution. 

The  resulting  radiation  will  be  the  extensively  studied  equilibrium 
air  radiation.  Hiis  equilibrium  air  radiation,  in  the  temperature 
range  of  interest,  has  a  very  steep  temperature  dependence,  so  that 
the  radiation  intensity  will  be  strongly  coupled  to  the  wake  temperature. 
Because  of  this,  the  decay  of  radiant  intensity  in  the  wake  will  be 
directly  related  to  the  temperature  decay  in  wake,  which  in  turn  is 
directly  related  to  the  decay  of  wake  enthalpy.  This  means  that  the 
wake  radiation  will  be  directly  related  to  the  overall  aerodynamic 
properties  of  the  wake. 

In  the  case  of  a  nonequilibrium  wake  there  is  very  little  atomic 
recombination  during  the  expansion  behind  the  body.  The  energy  of 
molecular  dissociation  is  not  liberated  during  this  expansion,  so  that 
the  temperature  drop  is  much  greater  than  for  the  equilibrium  case;  by 
the  time  the  pressure  has  decayed  to  the  ambient  value,  the  temperature 
has  decreased  to  a  value  of  the  order  of  1000°K.  Thus,  the  nonequilibrium 
wake  will  be  colder  than  the  equilibrium  wake.  However,  it  will  also 
be  much  longer  than  the  equilibrium  wake;  the  length  of  the  nonequilibritam 
wake  will  be  determined,  not  by  the  time  required  for  the  temperature 
or  wake  enthalpy  to  decay,  but  by  the  time  required  for  the  atomic 
species  to  recombine.  (At  high  altitudes,  200,000  ft  for  example,  this 
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recombination  time  can  be  several  seconds.)  Because  of  the  low 
temperature,  collisional  excitation  will  be  of  negligible  importance 
as  a  source  of  radiation.  However,  the  nonequilibrium  wake  is  a  highly 
reactive  mixture  of  atoms  and  molecules,  and  there  are  numerous 
chemiluminescent  reactions  which  can  be  a  significant  source  of 
radiation.  Since  these  reactions  are  exothermic,  their  rate  constants 
are  relatively  insensitive  to  temperature.  Because  of  this,  the  intensity 
of  the  wake  radiation  will  also  be  relatively  insensitive  to  the  wake 
temperature.  However,  since  the  radiation  is  a  chemical  phenomenon,  it 
will  be  greatly  affected  by  changes  in  the  chemical  composition  of  the 
wake.  This  chemical  composition  will  in  turn  be  affected  by  the  various 
chemical  reactions  that  occur  in  the  wake,  and  by  the  extent  to  which 
the  outer,  molecular-rich  inviscid  flow  is  mixed  into  the  inner,  atomic 
part  of  the  wake.  Thus,  in  the  nonequilibrium  case  the  wake  radiation 
will  be  directly  affected  by  the  mixing  processes  in  the  wake,  but  it 
will  not  be  affected  significantly  by  the  temperature  decay  in  the 
wake  or  the  enthalpy  decay. 

The  above,  very  brief  summary  of  the  properties  of  the  radiation 
from  nonequilibrium  wakes  will  be  expanded  in  the  following  sections. 
Section  II  presents  a  discussion  of  the  various  chemiliminescent 
reactions  of  importance  in  the  wake,  and  Section  III  discusses  the 
various  wake  phenomena  important  in  determining  the  amount  of 
radiation  produced  by  the  chemiluminescent  reactions.  In  Section  IV 
numerical  estimates  are  presented  of  the  amount  of  radiation  emitted 
from  the  wake,  and  Section  V  is  a  summary  and  discussion  of  the 


results. 
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II .  CHEMILUMINESCENT  REACTIONS 

As  indicated  in  the  Introduction,  the  wake  of  a  blunt  body  at 
altitudes  above  about  50,000  ft  will  be  relatively  cool,  of  the  order 
of  1000°K,  but  will  contain  large  numbers  of  unrecoiabined  atoms 
(N  and  0).  As  we  shall  see  in  the  next  section,  there  will  also  be 
significant  quantities  of  IK)  present.  For  a  mixture  of  this  temperature 
and  chemical  composition,  the  most  important  source  of  radiation  will 
be  chemiluminescence.  The  most  significant  chemiluminescent  reactions 
will  be 


N  +  N  + 

(M) 

-♦  N2  +  (M)  +  hv 

(1) 

N  +  0  + 

(M) 

>•  NO  +  (M)  +  hu 

(2) 

NO  +  0  + 

(M) 

-  NO2+  (M)  +  hu 

(3) 

where  a  third  body  (M)  may  or  may  not  participate  In  the  reaction. 

There  now  exists  a  great  deal  of  experimental  information  regarding 
these  reactions.  In  this  section  we  summarize  those  details  concerning 
these  reactions  of  importance  for  the  problem  of  wake  radiation. 

A.  LEWIS -RAYLEIGH  NITROGEN  AFTERGLOW 

Reaction  (1)  is  the  well  known  Lewis -Rayleigh  nitrogen  afterglow, 

the  spectrum  of  which  consists  principally  of  the  first  positive  bands 

of  N2.  Recent  kinetic  studies  of  this  reaction  show  that  for  pressures 

-3  -2 

between  1  and  10  mm  Hg  (10  to  10  atmospheres)  the  overall  reaction 

(2)* 

follows  a  two-body  rate  law  with  a  rate  constant 
* 

Young  and  Sharpless  state  that  this  value,  and  the  other 
ehemi luminescent  rate  constants  given  in  this  section,  are  probably 

accurate  to  a  factor  of  2. 


=  3.1  X  10  ^^cm^/sec 


For  pressures  outside  this  range  it  is  unclear  at  the  present  time 

as  to  whether  the  reaction  follows  a  two-body  or  three-body  rate 
* 

law.  However,  to  assume  the  two-body  law  holds  is  as  reasonable  as 

(2) 

any  other  assumption. 

The  spectrum  of  the  Lewis-Rayleigh  afterglow  consists  primarily 

of  the  N2  first  positive  bands,  with  additional  weak  emission  in  the 

Y  bands  in  the  visible  and  infrared,  and  in  the  Ljnnan-Birge-Hopfield 

system  in  the  vacuiam  ultraviolet.  The  strongest  bands  of  the  first 

(3)** 

positive  system  observed  in  the  afterglow  are: 


^v  =  5  sequence 
Av  =  4  sequence 


5370  A  to  5440  A 
5750  A  to  5870  A 


(11,8) 

band 

(6,3) 

band 

(4,2) 

band 

(3,1) 

band 

(2,0) 

band 

(2,1) 

band 

(1,0) 

band 

(0,0) 

band 

6250  A 
o 

6620  A 
7500  A 
7630  A 
7750  A 
8720  A 
8910  ^ 
10420  A 


It  should  be  pointed  out  that  the  existence  of  a  two-body  rate 
law  over  a  limited  region  of  pressure  does  not  necessarily  imply  a 
purely  two-body  reaction.  Many  examples  are  known  of  complicated  chain 
reactions  which,  over  a  limited  pressure  region,  obey  a  two-body  rate 
law.  It  is  very  likely  that  the  chemiltjrainescent  reactions  discussed 
in  this  section  are  of  that  type. 

"Me 

The  relative  intensities  of  the  various  first  positive  bands 
excited  by  reaction  (1)  are  not  the  same  as  for  the  first  positive 
emission  from  equilibrium  air. 


Reaction  (2)  produces  the  blue  nitric  oxide  afterglow,  consisting 

of  emission  in  the  and  6  bands  of  NO.  The  reaction  mechanism  is 

complicated,  and  results  in  an  overall  rate  constant  that  is  partly 

(4) 

two-body  and  partly  three-body.^  TSie  NO  3  system  is  entirely  three - 

(2) 

body,  with  a  rate  constant 


Kp  =  2,4  X  10  ^^cm^/sec 


Hie  NO  Y  system  is  part  two-body  and  part  three-body.  The  rate 

^  ^  (2) 

constants  are 


K  = 

Y 


- 18  8 

8.2  X  10  cm  /sec  (two-body) 


2.4  X  10  ^^cm^/sec  (three-body) 


Hie  M)  6  system  is  entirely  two-body,  with  a  rate  constant 


Kg  =»  1.1  X  10  ^^cm^/sec 


These  values  result  in  a  total  rate  constant  for  reaction  (2)  which 

(2 1 

is  part  two -body  and  part  three-body,  with  the  values^  ' 


-17  3 

1.9  X  10  cm  /sec  (two -body) 


-34  6 

2.6  X  10  cm  /sec  (three-body) 


The  relative  intensity  of  the  emission  in  the  various  vibration- 


rotation  bands  of  the  ^,y>  and  6  systems  has  been  measured,  and  is  shown 
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in  Tables  1,2,  and  3*^^^  With  respect  to  ground-based  or  airplane- 

based  observations  of  this  radiation,  it  should  be  noted  that  all  of 

o 

the  Y  S  radiation  lies  at  wavelengths  below  3000  A,  where  it 

will  be  absorbed  by  the  ozone.  However,  approximately  36  per  cent  of  the 

o 

radiation  from  the  P  system  lies  above  3000  A.  Using  this  number,  one 

may  obtain  an  effective  rate  constant  for  reaction  (2),  for  radiation 
o 

above  3000  A,  of 

=  0,9  X  10  ^^cm^/sec 


In  addition  to  the  P,y,  and  6  bands  of  NO,  the  recombination  of 

atomic  nitrogen  and  oxygen  is  known  to  produce  emission  in  the  NO 

o  o 

Ogawa  bands,  in  the  region  from  7800  A  to  9700  A,  and  in  the 

22+  o  (4) 

N0(G  TT  -  A  S  )  band  system,  in  the  region  around  12200  A.^  '  Quantitative 

information  conceraing  the  rate  constants  for  these  two  band  systems 

is  unavailable  at  the  present  time.  However,  based  on  the  presently 

2  2  + 

accepted  mechanism  for  the  NO  afterglow  radiation,  the  C  tt  -  A  S 

system  should  have  an  intensity  of  the  same  order  of  magnitude  as  the 

(4) 

two  body  part  of  the  y  system.' 

C.  NO2  CONTINUUM 

Reaction  (3)  produces  the  so-called  NO2  "continuum."  This 
emission  is  not  really  a  continuum,  but  has  a  fine  grained  banded 
structure,  due  to  the  complicated  chain  mechanism  producing  the  overall 
reaction, The  emission  spectrum  begins  at  3875  X,  peaks  at  about 

o  o  (6) 

6500  A,  and  extends  to  at  least  14000  A.  The  reaction  mechanism 
results  in  an  overall  rate  constant  that  is  two-body,  at  least  for 
pressures  above  about  0.1  mm  Hg  (10  ^  atmospheres). 


The  value 
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(Data 

Vibration  band 
6-1 

3- 1 

4- 2 
3-2 
2-2 
3-3 

2- 3 

3- 4 
2-4 

1- 4 

2- 5 
0-4 
1-5 
0-5 
1-6 
0-6 
0-7 
0-8 
0-9 
0-10 
0-11 


Table  1 
NO  3  SYSTEM 
taken  from  Ref.  2) 

o 

\(A) 

2018 

2144 

2187 

2232 

2288 

2327 

2382 

2428 

2488 

2552 

2602 

2623 

2676 

2751 

2806 

2889 

3039 

3203 

3381 

3578 

3795 


Relative  intensity 
0.62 
1.20 
0.65 
1.52 
0.98 
1.87 
0.88 
1.52 
1.44 
0.89 
2.54 
1.48 
0.56 
1.90 
0.53 
2.32 
3.60 
3.16 
2.34 
1.65 
0.94 


0-12 


4035 


0.40 
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Table  2 
NO  y  SYSTEM 

(Data  taken  frcsn  Ref.  2) 


o 

Vibration  band  X(A)  Relative  intensity 


3-0  1959  1.03 
3-1  2036  0.74 
3-2  2115  0.83 
0-0  2266  2.49 
0-1  2367  3.12 
0-2  2475  2.06 
0-3  2592  1.36 
0-4  2717  0.70 


0-5 


2855 


0.24 
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Table  3 
NO  8  SYSTEM 

(Data  taken  from  Ref.  2) 


Ltion  band 

\(A) 

Relative  i: 

0-0 

1914 

1.15 

0-1 

1985 

1.50 

CM 

1 

O 

2061 

2.53 

0-3 

2141 

2.44 

0-4 

2227 

1.59 

0-5 

2318 

0.68 

of  the  rate  constant,  for  emission  in  the  wavelength  region  from 
3875  A  to  14000  I  is^^^ 

K„  =  6.4  X  10  ^^cm^/sec 

O 


D.  EFFECT  OF  TEMPESATURE  ON  GHEMILI3MIHE SCENT  REACTIONS 

TSie  values  of  the  chani luminescent  rate  constants  given  in  the 
above  discussion  have  all  been  the  result  of  room  temperature  measure¬ 
ments.  The  temperatures  of  interest  in  the  wake  will  be  of  the  order 
of  1000°K  to  2000°K,  however,  so  that  the  question  arises  of  the 
temperature  dependence  of  the  rate  constants. 

At  the  present  time  there  is  no  clear  experimental  information 
concerning  this  temperature  dependence.  However,  there  does  exist  a 

reasonable  understanding  of  the  detailed  mechanisms  for  these  three 

(3-5) 

chemiluminescent  reactions,'’  and  it  is  well  established  that  for 

each  of  these  reactions  the  activation  energy  is  zero*^  This  means 

that  the  rate  constants  will  not  be  exponentially  dependent  on  temperature; 

if  there  is  a  temperature  dependence,  it  will  be  slight. 

Attempts  have  been  made  experimentally  to  determine  the 

(7  8 ) 

temperature  dependence  of  the  HO  +  0  reaction.'  ’  '  These  attempts 

were  based  on  high  temperature  shock  tube  measurements  of  the  NO^ 

"continuum,"  together  with  the  assianption  that  the  high  temperature 

radiation  and  the  low  temperature  afterglow  were  both  produced  by  the 

same  radiative  reccmbination  mechanism.  It  has  recently  been  suggested 

that  this  assumption  is  not  correct,  and  that  the  temperature  dependence 

(9) 

deduced  frcaa  it  is  not  justified.'  ^ 
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In  view  of  the  uncertain  situation  concerning  the  temperature 
depend^ce  of  the  chemiluminescent  rate  constants,  the  only  well- 
established  fact  being  that  the  dependence  on  temperature  will  be  slight, 
and  in  view  of  the  fact  that  the  room  temperature  rate  constants  are 
only  known  to  within  a  factor  of  two,  the  most  reasonable  course 
at  the  present  time  seems  to  be  to  use  the  room  temperature  values, 
without  attanptlng  a  dubious  extrapolation  to  higher  temperature.  The 
error  made  in  doing  this  will  most  probably  be  no  greater  than  the 
error  inherent  in  our  knowledge  of  the  room  temperature  values. 

Besides  any  effect  on  the  rate  constants,  the  temperature  should 
also  affect  the  spectral  distribution  of  the  various  afterglow 
reactions.  No  information  is  available  on  this  effect,  but  it  may 
turn  out  to  be  more  important  than  the  rate  constant  effect. 

E.  OTHER  CHEMILllMINESC^NT  BEACTIQN5 

In  addition  to  the  three  reactions  discussed  in  detail 
above,  a  number  of  other  chemiluminescent  reactions  will  occur 
in  the  wake.  The  most  extensively  studied  of  these  is  the  oxygen 
recombination  reaction 

0  +  0  +  (M)  ©2  +  (M)  +  hu 

which  produces  the  atmospheric  and  Herzberg  systaas  of  0^.  This 

reaction  will,  for  two  reasons,  be  of  much  less  importance  than  the 

previous  three  as  a  source  of  wake  radiation:  (1)  The  rate  constants 

are  much  smaller;  for  the  atmospheric  bands  a  three-body  rate  constant 

of  the  order  of  10  ^^cm^/sec; for  the  Herzberg  bands  a  two-body 

-21  3  (2) 

rate  constant  of  the  order  of  10  cm  /sec',  '  (2)  The  radiative 

transitions  involved  are  partially  forbidden  with  relatively  long 


u 


radiative  lifetimes,  and  appear  to  be  strongly  quenched  by  collisional 

(2) 

de-excitation. 

o  o 

The  red  (6300  A)  and  green  (5577  A)  lines  of  atomic  oxygen  are 

also  produced  by  chemilmaine scent  reactions.  Hie  red  line,  with  a 

-2  -1 

radiative  rate  of  about  10  sec  ,  will  certainly  be  collisionally 

de-excited  at  reentry  altitudes.  Hie  situation  with  respect  to  the 

green  line  is  less  clear.  Also,  the  absolute  rate  of  chemiluminescent 

(2) 

production  of  the  green  line  is  uncertain  at  the  present  time. 

In  addition  to  the  above  reactions,  which  produce  radiation  in 
the  UV,  visible,  and  near  IR,  numerous  chemiluminescent  reactions 
produce  radiation  in  the  middle  IR. The  most  important 
spectral  feature  in  this  region,  for  a  pure  air  wake,  will  probably  be 
the  fundamental  vibration- rotation  band  of  NO  at  5.3  microns,  and 
the  most  important  reaction  producing  this  band  will  probably  be 

N  +  O2  ^  NO*  +  0 

where  the  superscript  *  denotes  a  vibrationally  excited  molecule  in  its 
ground  electronic  state.  The  amount  of  vibrational  excitation  produced 
by  this  reaction  has  not  yet  been  measured.  However,  the  similar 
reaction 

N  +  NO  -»  N*  +  0 

is  known  to  produce  large  amounts  of  vibrationally  excited 
Because  of  the  low  excitation  energy  of  the  NO  vibration-rotation  band, 
about  0.2  ev,  and  because  of  the  long  radiative  lifetime,  of  the  order 
of  10  ^sec,  collisional  excitation  will  also  be  an  Important  source 
of  NO*,  even  at  altitudes  of  150,000  to  200,000  ft. 
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F.  FACTORS  AFFECTOG  CTOIILIJHINSSOTT  RADIATION 

As  we  have  discussed  above,  the  various  chemi luminescent  reactions 
expected  to  be  of  importance  as  sources  of  wake  radiation  are  not  very 
temperature  sensitive.  'Riis  is  in  contrast  to  the  radiation  from 
equilibrium  air,  which  is  exponentially  dependent  on  temperature.  The 
amount  of  chemiluminescent  radiation  produced  will,  however,  be  very 
dependent  on  the  species  concentrations  in  the  wake.  It  will  depend 
on  the  amount  of  N,0,  and  NO  present  in  the  wake,  and  it  will  depend 
on  the  length  of  time  that  these  species  are  present.  In  the  following 
section  we  will  discuss  the  various  wake  phencxaena  that  are  important 
in  determining  the  species  concentrations  in  the  wake,  and  therefore 
also  important  in  determining  the  amount  of  wake  radiation. 


is 

III>  WAKE  ra£EK)M^LQGY 

As  indieated  in  the  previous  section,  the  amotmt  of  radiation 

from  a  nonequilibrium  wake  will  depend  primarily  on  the  species 

concentrations  in  the  wake,  in  particular  on  the  concentrations  of 

N,0,  and  NO,  since  these  are  the  reactants  involved  in  the  most 

important  chemiluminescent  reactions.  The  value  of  these  species 

concentrations  will  in  turn  depend  on  a  number  of  flow  field  phencHaena: 

the  conditions  established  in  the  stagnation  region,  the  expansion  of 

the  flow  around  and  immediately  behind  the  body,  the  chemical  reactions 

•k 

occuring  in  the  far  wake,  and  the  effect  of  turbulent  mixing  in  the 
wake.  In  this  section  we  will  discuss  the  effect  of  these  various 
phenomena  on  the  species  concentrations. 

A.  STAGNATION  AND  EXPANSION  REGION 

Upon  crossing  the  bow  shock  of  a  blunt  body,  the  air  flow  is 
rapidly  compressed  and  heated  to  a  very  high  temperature.  Molecular 
dissociation  occurs,  and  the  chemical  composition  of  the  air  begins  to 
approach  equilibrlxan.  Whether  or  not  chonical  equilibrium  is  achieved 
in  the  stagnation  region  at  the  front  of  the  body  depends  on  a  number 
of  factors:  the  free  stream  velocity,  the  ambient  density,  and  the 
nose  radius  of  the  body.  For  typical  ICSM  reentry  velocities,  of  the 
order  of  20,000  ft/sec,  and  for  reentry  bodies  of  typical  size,  with  nose 
radii  of  the  order  of  1  ft,  it  has  been  found  that  the  stagnation 

(12) 

region  will  be  in  equilibrium  for  altitudes  up  to  about  200,000  ft.' 
Assuming  that  equilibrium  is  attained,  the  mole  fraction  species  con¬ 
centrations  in  the  stagnation  region  will  be  as  shown  in  Figs.  1-3. 

By  "far  wake"  we  will  mean  that  part  of  the  wake  at  ambient 
pressure. 
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Altitude  (ft) 

Fig.  1 — Atomic  nitrogen  concentration  in  stagnation  region 
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Altitude  (ft) 


Fig. 2— Molecular  nitrogen  concentration  in  stagnation  region 
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(13) 

The  data  in  Figs.  1-3  were  obtained  from  the  Avco  gas  dynamic  charts, 
together  with  Gilmore's  equilibrium  air  tables. 

As  the  airflow  leaves  the  stagnation  region,  it  begins  to 
expand  around  and  behind  the  body.  This  expansion  continues  until  the 
pressure  has  decayed  to  the  free  stream  value,  which  occurs  at  a  down¬ 
stream  distance  X/R  f^bTM.S,  where  R  is  the  nose  radius  and  M 
the  free  stream  Mach  number. During  this  expansion  the  air 
cools  rapidly,  and  for  chemical  equilibrium  to  be  maintained,  atomic 
reccaabination  should  occur  in  this  expansion  region.  However,  for 
low  enough  ambient  density  (or  high  enough  altitude)  the  atomic 
recombination  processes  will  lag  behind  the  temperature  drop,  resulting 
in  a  nonequilibrium  chemical  composition. 

Recent  studies  of  the  cheaaically  reacting,  expanding  flow  behind 
blunt  bodies  show  that  this  nonequilibrium  effect  will  be  of  major 
importance  for  altitudes  above  about  50,000  ft.^^^  Earlier  studies 
have  shown  that  for . altitudes  above  100,000  ft,  the  mole  fraction 
concentrations  of  the  major  species  in  the  flow  (N,  N2 »  0)  freeze  at 
essentially  their  stagnation  point  values . However,  and  this 
will  be  of  importance  in  the  later  discussion,  those  species  present 
in  the  stagnation  region  only  in  small  amounts  (O25  M))  do  not  freeze 
during  the  expansion;  due  to  the  effect  of  certain  bimolecular 
exchange  reactions  to  be  discussed  later,  the  concentrations  of  O2 
and  NO  are  reduced  by  over  an  order  of  magnitude  below  the  values  that 
would  be  estimated  from  any  sudden  freezing  approximation.^  ^ 

* 

The  so-called  P  P  point. 
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Since  those  species  carrying  the  bulk  of  the  chemical  enthalpy 
(N  and  0)  freeze  during  the  expansion,  the  temperature  drop 
associated  with  the  expansion  is  much  greater  than  for  an  equilibrium 
x^ake,  B€C:aus€  the  chmaical  composition  is  essentially  constant  during 
the  expansion,  the  expansion  process  will  be  approximately  isentropic, 
and  the  value  of  temperature  and  density  at  the  end  of  the  expansion 
region  (the  P  -♦  point)  can  be  estimated  from  the  usual  isentropic 
relations  T  oc  P^*^  p  a  P^^^,  where  y  is  now  the  effective 

(19) 

specific  heat  ratio  appropriate  to  the  frozen  chemical  ccmposltlon. 

Using  this  method,  together  with  the  mole  fraction  values  from 
Figs.  1-3,  one  obtains  the  value  of  temperature  and  density  at  the 
P  -»  P  point  shown  in  Figs.  4  and  5.  The  concentrations  of  the 
frozen  species  N,  H2,  and  0  at  the  P  -♦  P^  point  can  then  be  obtained; 
these  are  shown  in  Figs.  6-8.  These  results,  although  obtained  by 
assuming  a  simple  isentropic  expansion,  are  in  substantial  agreeiBent 
with  results  of  more  detailed  numerical  calculations . 

Several  comments  are  now  in  order  concerning  these  results.  First 
of  all,  the  values  of  traaperature  shown  in  Fig.  4  bear  out  our  previous 
remarks;  the  far  wake  will  he  relatively  cool,  with  temperatures  initially 
of  the  order  of  1000°K.  For  temperatures  of  this  magnitude,  as  we  have 
discussed  previously,  the  important  chemical  reactions  will  be  exothermic 
ones  which  are  relatively  insensitive  to  temperature.  In  addition. 

Fig.  3  shows  that  -Khen  the  flow  has  reached  the  P  -♦  P^  point,  the  density 
will  be  within  a  factor  of  approximately  four  to  seven  of  the  free 
stream  value.  Downstream  of  the  P  point,  the  temperature  and 

density  will  asymptotically  approach  the  free  stream  values.  Since 

Vibration  equilibrium  was  assumed  during  the  es^pansion  process, 
since  this  has  been  shown  to  be  a  good  approximation,'^®) 
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Fig.  6— Atomic 
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Fig8~Atomic  oxygen  concentration  at  P-»^Pco  point 


they  are  not  too  far  removed  frcrai  these  values  to  begin  with,  the 
temperature  and  density  variations  in  the  far  wake  will  be  relatively 
small,  much  smaller  than  for  the  case  of  an  equilibrium  wake* 

The  results  shown  in  Figs.  4-8  are  the  values  on  the  axis.  The 

off-axis  behavior  of  these  quantities  will  also  be  of  importance. 

Itetailed  numerical  studies have  shown  that,  for  the  special  cases 

considered,  the  concentrations  of  the  atomic  species  N  and  0  will 

r^aain  relatively  constant  for  some  distance  off-axis,  and  then  decrease 

rapidly  to  zero.  The  temperature  will  experience  scaae  off-axis  variation, 

but  will  remain  within  the  range  of  1000°K  to  2000°K  for  about  the  same 

distance  off-axis,  and  then  will  also  drop  rapidly. These  detailed 

studies  suggest  that  a  reasonable  model  of  the  nonequilibrium  wake 

may  be  obtained  by  assuming  a  cylindrical  region  surrounding  the  axis 

over  which  the  axis  values  of  species  concentration  (N,  N2  and  0)  and 

temperature  hold.  The  radial  extent  of  this  region  varies  with  down- 

(1) 

stream  position;  at  the  P  -*  P  point  it  is  of  the  order  of  3  R  , 

”  n 

Hiere  is  one  exception  to  the  above  statement  concerning  the 
off-axis  behavior  of  the  species  concentrations.  13ie  NO  concentration 
is  very  low  on  the  axis,  but  there  is  a  thin  shell  of  relatively  high 
NO  concentration  about  3  off-axis. The  significance  of  this  NO 
shell  will  be  discussed  later. 

The  above  discussion  pertains  to  the  values  attained  by  the  various 
wake  parameters  at  the  end  of  the  expansion  region.  Downstream  of  this 
point  the  wake  will  continue  to  cool,  but  since  the  temperature  is  so 
low  to  begin  with,  the  cooling  will  not  greatly  affect  the  chemical 
ccaaposltion.  The  chemical  ccwipositi on  will  be  affected,  however,  by 
the  various  chmical  reactions  that  occur  in  the  wake,  and  by  the 


mixing  processes  that  occur.  We  will  now  discuss  these  effects. 


B-  CHEMICAL  KINETICS  Of  FAR  WAKE 

As  shown  above,  the  initial  concentrations  of  N,  and  0  in  the 
far  wake  (i.e.,  the  concentrations  at  the  P  -*  point)  are  determined 
by  the  processes  occurring  during  the  expansion.  The  duration  of  the 
N  and  0  in  the  wake,  either  in  terms  of  lifetime  or  downstream  distance, 
will  be  determined  by  the  chi^ical  reactions  that  occur  in  the  wake. 
These  reactions  will  also  determine  the  amount  of  0^  and  liO  present 
in  the  wake. 

TSxe  detailed  nature  of  the  ch^ical  kinetics  in  the  wake  will  be 
strongly  affected  by  the  amount  of  mixing  that  occurs  between  the 
inner  wake,  which  is  mostly  atomic,  and  the  outer,  molecular  part  of 
the  wake.  For  the  moment  we  will  neglect  this  mixing;  later  in  this 
section  we  will  discuss  its  effect. 

To  begin  with,  we  consider  the  0^  kinetics.  In  the  absence 
of  mixing,  the  in  the  inner  part  of  the  wake  will  be  produced  by 
the  reaction 

0  +  0  +  M-*  O2  +  M  (4) 


Because  of  the  chemical  composition  of  the  wake,  i.e.,  large  quantities 
of  atomic  oxygen  present,  and  because  of  the  enhanced  catalytic  efficiency 
of  0  relative  to  N  or  N2»  0  will  be  the  most  effective  catalyst  for  this 

O 

reaction,  i.e.,  M  s  0,  so  that  the  rate  of  reaction  (4)  will  be  F,  =  K.N;:, 

4  4  0 

where  is  the  concentration  of  0,  and  the  rate  constant 

=  2,1  X  lO’^V^^^cm^/sec 


where  T  is  in  degrees  Kelvin,  The  0^  will  be  consumed  primarily  by 
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N  +  0^  M)  +  0 


with  a  rate  ,  and  a  rate  constant 


(5) 


(21) 


1.  O  O  a.  -3540/T  3, 

Kc  =  2.2  X  10  T  e  an  /sec 

o 


A  steady  state  is  very  qtiickly  reached  between  reaction  (4)  producing 
O2  and  reaction  (5)  constanlng  it.  The  steady-state  concentration  of 


©2  will  be 


N, 


Vo 


This  steady  state  is  achieved  in  a  time  ~  l/K^K^^,  which,  for  the 

temperature  and  nitrogen  atom  density  initially  present  in  the  wake, 

-3 

is  of  the  order  of  10  seconds. 

The  IK)  in  the  inner  wake,  in  the  absence  of  mixing,  will  be 
produced  by  reaction  (5),  and  by 


N+O+M  NO+M 


(6) 


which  has  a  rate  F^  =  and  a  rate  constant 


(23) 


K,  =  2.8  X  10'^®  T'^^^cm^/sec 
o 


for  all  catalysts  except  KO.  (Since  the  NO  concentration  will  be  very 
low,  it  will  not  be  an  important  catalyst.)  The  NO  will  be  consimaed  by 


N  +  NO  N2  +  0 


(7) 


which  has  a  rate  F^  ==  and  a  rate  constant 


(21,22) 


3© 


=  2.5  X  10  cm^/sec 


13iese  reactions  rapidly  produce  an  NO  steady-state  concentrar€ion 
given  by 

^0  K, 


This  steady-state  concentration  is  attained  in  a  time  ^  l/K^N^^, 
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which  for  the  cases  of  interest  here  is  of  the  order  of  10  seconds, 
If  the  ©2  steady  state  is  also  attained,  we  will  have 


N, 


NO 


■^0  •^Wo*S. 


Because  the  0_  and  NO  chemical  relaxation  times,  and 

2  O2  NO 

are  so  short,  the  O2  NO  concentrations  are  not  frozen  in  the  wake, 
but  remain  in  a  quasi -equilibrium  with  respect  to  the  N  and  0 
concentrations.  These  short  relaxation  times  are  due  to  the  fast 
bimolecular  exchange  reactions,  reactions  (5)  and  (7).  The  importance 
of  these  reactions  depends  on  the  presence  of  atomic  nitrogen,  and 
they  will  maintain  the  quasi-equilibrium  of  O2  and  NO  as  long  as  N  is 
present  in  significant  quantities.  These  same  bimolecular  exchange 
reactions  are  responsible  for  the  fact  mentioned  previously,  that  the 


O2  and  NO  do  not  freeze  in  the  expansion  region. 

We  can  now  consider  the  N  and  0  kinetics,  which  will  involve  much 
longer  characteristic  times  than  the  ©2  or  NO  kinetics.  As  the  N  and 
0  concentrations  change,  the  O2  and  NO  concentrations  will  quickly 
adjust  to  follow  the  change,  as  long  as  N  is  present. 
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Let  us  now  discuss  the  N  kinetics.  Atomic  nitrogen  is  not 
being  produced  by  any  reactions.  It  is  being  consinaed  by  reactions  (5), 
(6),  and  (7),  and  by 

N  +  N  +  M  ^  H2  +  M  (8) 

As  long  as  significant  quantities  of  atcanic  nitrogen  are  present,  N 
will  be  the  most  important  catalyst  for  this  reaction,  so  that 

the  order  of  1000°K,  is  rather  uncertain  at  the  present  time.  Wray^^^^ 
gives 

„  ,  ,  -3l2  6, 

Kq  =  6.4  X  10  T  cm  /sec 

O 

which  for  T  1000°K  gives  Kg  =  2  x  10  ^^cm^/sec.  Bortner^^^^  gives 

„  o  1 n~29  _-l  6 , 

Kg  =  2  X  10  T  cm  /sec 

r\  —  *^9  fi 

which  for  T  1000  K  gives  Kq  =  2  x  10  cm  /sec.  TSiere  is  thus  a 

O 

factor  of  ten  uncertainty.  Both  of  these  values  are  obtained  from 
high  temperature  shock  tube  experiments,  the  results  of  which  are 
extrapolated  to  lovrer  tanperature.  There  is  only  one  direct  experimental 
measurement  available  in  the  tanperature  region  of  interest  here,  that 
of  Back  et  al.,  at  700°K.^^^^  They  found  that  was  of  the  order  of 

O 

"3X  6 

(1-3)  X  10  cm  /sec.  Based  on  this  measurement,  we  will  here  use 
Wray's  value,  although  the  situation  is  still  somewhat  uncertain. 

For  the  range  of  values  of  wake  parameters  of  interest  here,  all 
four  of  the  reactions  consuming  N  can  be  of  importance,  so  that  the 
atcmalc  nitrogen  recombination  time  will  be  given  by  ~  N^^/ (r^+rg+r2+2rg) , 
or 


KgU^.  The  rate  constant  for  this  reaction,  for  temperatures  of 
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_ 1 _ 

Vo%  ^  ^>^*4 

This  is  approximately  the  length  of  time  that  atomic  nitrogen  will 
persist  in  the  wake.  While  it  is  present,  there  will  be  negligible 
oxygen  recombination,  because  reactions  (5)  and  (7)  quickly  liberate 
any  0  removed  by  reactions  (4)  and  (6).  Once  the  N  is  removed,  reaction 
(4)  will  begin  to  consume  0,  and  the  O2  concentration  will  rise,  "flie 
atomic  oxygen  recombination  time  will  be  approximately 


T 


0 


1 


As  later  numerical  results  will  show,  for  all  cases  of  interest 

we  will  have  T„  «  T^,  so  that  there  will  be  a  period  of  time  during 
N  0 

which  the  wake  will  consist  mostly  of  H2  and  0,  i.e,,  very  little  N. 

During  this  intermediate  period  (after  N  is  gone,  before  0  is  gone) 

the  NO  concentration  will  be  the  same  as  before,  since  both  the  NO 

production  and  destruction  reactions  depend  on  Nj^  in  the  same  manner. 

The  only  difference  will  be  that  the  NO  relaxation  time,  =  I/K^N^^, 

will  become  much  larger,  i.e.,  the  NO  concentration  will  "freeze"  at 

* 

the  value  established  virile  N  was  present.  This  frozen  NO  concentration 
will  persist  for  a  time  much  greater  than  the  oxygen  recombination  time, 
and  the  ultimate  fate  of  the  NO  need  not  concern  us  here. 

The  reaction  NO  +0  +  M  NO2  +  M  will  not  affect  the  IK) 
concentration,  since  it  will  be  inBaediately  followed  by  NO2  +  0  NO  +  ©2 

which  recreates  the  M©  molecule. For  the  N©  concentrations  achieved 
in  the  wake,  these  two  reactions  will  also  have  a  negligible  effect  on 
the  oxygen  recombination  rate. 
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For  the  problem  of  determining  the  amount  of  wake  radiation  the 

important  quantities  determined  by  these  various  chemiical  reactions  are 

the  concentration  of  IK)  in  the  wake,  the  duration  of  N  and  0 

in  the  wake,  and  Tq.  In  Figs.  9  and  10  we  plot  these  quantities, 

using  the  axis  values  of  the  concentrations  of  N,  N2,  and  0,  l.e,, 

the  values  shorn  in  Figs.  6-8,  since  these  axis  values  are  appropriate 

for  the  inner  wake.  In  determining  the  rate  constants,  we  used  a 
o 

temperature  of  1000  K,  which  is  roughly  appropriate  for  the  cases  of 

•ft 

interest  here.  The  values  shown  in  Figs.  9  and  10  apply  to  the  inner 
wake,  i.e.,  the  high -enthalpy ,  atomic-rich  part  of  the  wake.  As  we 
said  before,  they  are  correct  in  the  absence  of  mixing  between  the 
inner  and  outer  wake. 

C.  EFFECT  OF  TUKBULEUT  MIXING 

We  must  now  consider  the  effect  of  mixing  between  the  inner  and 
outer  wake  on  the  chemical  kinetics  and  species  concentrations  in  the 
wake.  One  effect  of  this  mixing  will  be  to  cool  the  inner  wake.  This 
effect  has  been  extensively  studied  for  botli  laminar  and  turbulent’ 

(Of. 

mixing,''  ’  ^  However,  since  the  radiation-producing  reactions  of 

interest  to  us  are  relatively  insensitive  to  temperature,  and  the 
temperature  is  low  to  begin  with,  this  cooling  effect  will  be  of  no 
great  significance  as  far  as  the  nonequilibri\im  radiation  is  concerned. 
What  will  be  Important  is  the  extent  to  vdilch  the  mixing  process  changes 
the  chemical  ccmaposition  of  the  wake,  or  more  specifically,  the  extent 
to  which  it  affects  the  NO  concentration  in  the  wake,  and  the  nitrogen 
and  oxygen  recombination  times. 

Since  the  rate  constants  are  not  very  sensitive  to  te»perature 
changes,  using  one  value  of  tmperature  does  not  introduce  a  significant 
error. 
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Fig.9— NO  concentration  in  inner  wake  (no  mixing) 


The  steady-state  N©  coaeeutxation  la  the  wake  is  established 
during  the  initial,  atomic -nitrogen -rich  part  of  the  wake.  ISie  chemical 


reactions  occurring  during  this  period  also  determine  the  nitrogen 

reccmbination  time.  The  chemical  kinetics  of  this  atcmaic  nitrogen- 

•k 

rich  region  are  strongly  affected  by  the  low  0^  concentration, 

which  limits  the  rate  at  which  reaction  (5)  produces  NO,  and  which 

also  helps  to  limit  the  rate  at  isdiich  reactions  (5)  and  (7)  remove  N. 

(1) 

As  pointed  out  by  Lin  and  Hayes,'  the  mixing  of  additional  0^  from 
the  outer,  molecular-rich  part  of  the  wake  into  the  inner,  atomic- 
nitrogen-rich  wake  can  have  a  profound  effect  on  the  wake  kinetics, 
increasing  the  NO  concentration,  and  decreasing  the  nitrogen  reccanbina- 
tion  time.  The  magnitude  of  the  effect  will  depend  on  the  rapidity  of 
the  mixing  process. 

If  the  wake  were  to  remain  laminar  down  to  low  altitudes,  so  that 

mixing  occurred  only  by  molecular  diffusion,  the  effect  of  mixing  on 

the  wake  kinetics  would  not  be  very  large.  This  is  because,  for  reentry 

bodies  of  typical  size  (R^  ^  1  ft),  the  characteristic  laminar  diffusion 

f  28^ 

time  is  ^  1  sec  for  altitudes  below  about  200,000  ft. ^  ^  For  typical 

ICSiN  velocities  (7  to  8  km/sec)  this  is  greater  than  the  nitrogen 

reccmblnation  time  (cf.  Fig.  10),  so  that  by  the  time  the  ©2  diffuses 

into  the  inner  ^ke  the  N  has  already  recombined. 

However,  the  wake  will  not  remain  laminar.  Recent  data  indicate 

that  for  blunt  bodies  the  entire  wake  downstream  of  the  neck  will  be 

4  (29) 

turbulent  for  freestream  Reynolds  numbers  above  about  5  x  10  . 


!Rie  0_  to  0  ratio  in  the  absence  of  mixing  is  which 

is  of  the  order  of  10“^^  to  10“^. 

Mixing  of  the  NO  shell  into  the  inner  wake  can  also  affect  the 
kinetics.  However,  nun^rlcal  studies''  '  have  shown  this  to  be  of  lesser 
importance  than  the  ©2  mixing,  and  we  will  not  consider  it  further  here. 
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For  a  1-ft  body  at  Mach  20  dils  corresponds  to  altitudes  up  to  about 

200,000  so  that  over  the  altitude  region  of  interest  here,  the 

entire  wa’ke  will  be  turbulent.  Turbulent  diffusion  processes  are 

(27) 

known  to  be  much  more  rapid  than  laminar  ones,  so  that  in  the 
presence  of  turbulence  it  is  possible  that  sufficient  0^  could  be 
mixed  into  the  inner  wake  to  affect  the  kinetics. 

Unfortunately,  at  the  present  tia^  there  are  large  uncertainties 
in  our  knowledge  of  the  turbulent  diffusion  process.  In  particular, 
it  is  not  known  how  rapidly  molecular  mixing  will  occur  within  the 
turbulent  region.  IHiis  means  that  we  cannot  specify  with  any  assurance 
the  rate  at  which  0^  from  the  outer  wake  will  be  microscopically  mixed 
into  the  inner  wake. 

In  view  of  this  uncertainty,  Lin  and  Hayes  have  proposed  two 

extreme  models  for  the  turbulent  mixing  process,  which  it  is  hoped 

(1) 

will  bracket  the  actual  physical  situation.'  ^  In  the  first  of 
these  models,  termed  "inviscid  random  convection,"  the  turbulent 
eddy  motion  is  assumed  to  mix  cold  lumps  of  air  from  the  outer  wake 
into  the  hotter  inner  wake,  so  that  the  turbulent  region  consists  of 
a  mixture  of  hot  and  cold  lumps.  In  this  inviscid  randcmi  convection 
model  the  molecular  diffusion  between  hot  and  cold  limaps  is 
assumed  negligible,  so  that  a  large-scale  mixing  occurs  between 
the  inner  and  outer  parts  of  the  wake,  but  no  microscopic  mixing. 

* 

As  an  example  of  this  uncertainty.  Ref.  30  gives  four  turbulent 
diffusivity  models,  all  of  which  are  plausible  on  the  basis  of  existing 
information,  but  which  produce  orders  of  magnitude  difference  in  the 
values  of  the  various  wake  parameters. 
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Since  there  is  no  miGroscopic  mixing,  the  choice  1  kinetics  of  the 
fluid  elements  will  he  the  same  as  for  the  non- turbulent  case,  so  that 
the  quantities  of  particular  interest  to  us,  the  NO  concentration  and 
the  nitrogen  and  oxygen  recombination  times,  will  be  unaffected  by 
the  turbulence,  i.e.,  they  will  still  have  the  values  shown  in  Figs.  9 
and  10. 

The  second  model  of  turbulence  proposed  by  Lin  and  Hayes,  the 
homogeneous  mixing  model,  is  intended  to  represent  an  upper  limit 
on  the  effect  of  turbulence.  In  this  model  the  entire  turbulent 
region  is  assumed  to  be  homogeneous,  and  microscopic  mixing  of  the 
new  fluid  crossing  the  turbulent  boundary  with  the  old  fluid  already 
within  the  boundary  is  assumed  to  occur  instantaneously.  The 
entrainnmnt  of  02-rich  air  by  the  turbulent  boundary  is  therefore 
equivalent,  in  this  model,  to  a  uniform  volume  production  rate  of  O2 
within  the  inner  wake. 

We  will  use  this  homogeneous  mixing  model  of  Lin  and  Hayes  to 
obtain  an  upper  limit  on  the  effect  of  turbulence  on  the  wake  parameters 
of  interest  here.  To  begin),  we  must  estimate  the  amount  of  ©2 
injected  into  the  inner,  turbulent  wake.  The  amount  of  crossing 
the  turbulent  entrainment  boundary  per  unit  length  of  wake  will  be 


given  by 


2Tr 

t 


dY 

f 

dX 


u 

e 


where  X  is  the  axial  distance,  Yj(X)  is  the  radius  of  the  turbulent 

front,  and  and  u  are  the  0_  concentration  and  stream  velocity 

in  the  inviscid  flow  just  outside  the  turbulent  core.  The  volume 

2 

of  the  turbulent  core  per  unit  length  of  weke  Is  tt  Y^,  so  that. 


assuming  instantaneous  homogeneous  mixing,  the  equivalent  volume 
production  rate  of  0^  is 


u 

e 


For  X/R  ^  10^,  where  R  is  the  nose  radius,  the  radius  of  the 
n  ’  n 

turbulent  front  is  approximately^^^ 


0.8 


1/3 


^iSiis  gives  dY^/dX  =  Y^/SX,  so  that 


2 


3X 


We  are  interested  in  the  entrainment  of  the  02~  rich  part  of  the 
inviscid  flow.  This  will  begin  to  occur  once  all  of  the  dissociated 
inviscid  flow  has  been  engulfed  by  the  turbulent  front.  According 

Cl) 

to  calculations  of  Lin  and  Hayes,'  ^  this  point  is  reached  somewhere 

between  100  and  600  nose  radii  downstream,  depending  on  the  body  size 

and  altitude.  The  02“ rich  part  of  the  inviscid  flow  will  be  cool, 

with  temperatures  of  the  order  of  1000°K  or  less,  and  will  be  at 

ambient  pressure,  so  that  ,  the  ambient  value  of  O2  concentre- 

*  2  2 

tlon.  The  flow  in  this  region  will  be  at  essentially  free  stream 
velocity,  i.e.,  u^  Finally,  for  bodies  of  about  1-ft  size 

we  are  interested  in  values  of  X  ^  10  ft.  Since  u  ~  10^  ft/sec,  we 
obtain  as  an  order  of  magnitude  estimate  of  the  ©2  injection  rate 
* 

Tnis  may  slightly  overestimate  the  0^  concentration. 
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10 


-1 

sec 


This  O2  injection  into  the  inner  wake  will  now  modify  the  wake 
kinetics.  It  will  provide  a  source  of  O2  in  addition  to  reaction  (4), 
so  that  the  steady  state  O2  concentration  (in  the  atomic  nitrogen  part 
of  the  wake)  will  now  be 


k,n:  +  P- 
4  0  0- 


KjN^ 


since  for  the  cases  of  interest  here,  the  turbulent  injection  will 
dominate  over  reaction  (4).  The  steady  state  NO  concentration  will 


be 


^0 


K  N_ 
5  O': 


and  the  nitrogen  recombination  time  will  become 


■nie  oxygen  recombination  time  will  be  left  essentially  unchanged  by 

the  turbulence;  the  atomic  oxygen  will  be  spread  out  over  a  larger 

region,  but  it  will  still  recombine  by  reaction  (4), 

In  Figs,  11  and  12  we  present  the  values  of  and  for  this 

instantaneous  homogeneous  mixing  model.  By  comparison  with  Figs.  9 

and  10,  we  see  that  the  presence  of  turbulence  has  increased  the  NO 

4 

concentration  in  the  inner  wake  by  a  factor  of  between  10  and  10  , 


Particles  /cm 


A1 


Velocity  (ft/sec) 

Fig.  11— NO  concentration  in  wake  (homogeneous  mixing) 


Velocity  (ft/sec) 

Fig.  12 — Nitrogen  recombination  time  (homogeneous  mixing) 
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depending  on  the  altitude  and  velocity.  Except  for  the  case  of  lowest 
altitude  (100,000  ft)  and  highest  velocity  (24,000  ft/sec),  the  turbulent 
mixing  has  also  decreased  the  nitrogen  recombination  time.  At  the 
higher  altitudes  (200,000  ft)  this  decrease  can  be  as  much  as  four 
orders  of  magnitude.  The  large  magnitude  of  this  effect  is  due  to 
the  fact  that,  compared  to  the  three-body  recombination  of  atomic 
oxygen,  the  turbulent  mixing  of  ©2  into  the  inner  wake  is  a  very 
efficient  source  of  02»  especially  at  the  higher  altitudes. 

The  values  of  T„  shown  in  Fig.  12  represent  the  length  of  time 
that  atomic  nitrogen  will  remain  in  the  wake  once  the  turbulent  front 
has  begun  to  engulf  the  undlssoclated  flow.  This  will  begin  to  occur 
a  few  hundred  nose  radii  downstream.  For  typical  ICBM  velocities  and 

a  nose  radius  of  the  order  of  1  ft,  this  will  correspond  to  a  time 

-2 

of  the  order  of  10  sec.  This  is  larger  than  the  values  shown  in 

Fig.  12,  so  that  for  the  case  of  homogeneous  mixing  the  total  length 

of  the  atomic  nitrogen  wake  will  be  determined  primarily  by  the  time 

required  for  the  turbulent  core  to  expand,  and  not  by 

As  mentioned  previously,  the  inviscid  random  convection  model 

and  the  homogeneous,  mixing  model  are  thought  to  represent  the  two 

possible  extremes  for  the  turbulent  mixing  in  the  wake,  with  the  actual 

physical  situation  lying  somewhere  in  between.  The  results  shown  in 

Figs.  9  and  10  for  no  mixing,  and  in  Figs.  11  and  12  for  homogeneous 

mixing  therefore  represent  the  possible  range  of  values  for 

and  T  .  Because  of  the  uncertainties  concerning  the  nature  of  the 
H 

turbulent  process,  it  is  not  possible  at  present  to  delimit  these 
quantities  with  any  greater  accuracy. 


l^ere  is  one  further  effect  of  turbulence  which  should  be  conoaented 
on.  As  pointed  out  by  Lin  and  Hayes, the  enhancement  of  the  exo¬ 
thermic  reactions  (5)  and  (7)  when  additional  0^  molecules  arc  entrained 
into  the  wake  from  the  external  flow  results  in  a  reheating  of  the 
wake.  iSieir  results  show  that  for  some  altitudes  the  magnitude  of  this 
reheating  may  be  as  large  as  1000°K.  This  effect  will  increase  the 
wake  temperature  above  the  value  used  here,  but  not  by  a  large  enough 
amount  to  affect  our  conclusions.  The  rate  constants  of  the  important 
chemical  reactions  will  not  be  altered  by  a  large  amount,  and  chemi¬ 
luminescent  reactions  will  still  be  the  principal  source  of  radiation. 


Having  discussed  the  important  chemiluminescent  reactions,  and 
the  important  chemical  and  aerodynamic  phenomena  affecting  the  chemical 
composition  of  the  wake,  we  are  now  in  a  position  to  calculate  the 
amount  of  radiation  fflaitted  by  the  nonequilibrium  wake  of  a  reentry 
vehicle.  In  carrying  out  this  calculation,  we  will  be  primarily 
concerned  with  obtaining  an  order  of  magnitude  estimate  of  the  radiant 
intensity,  determining  the  altitude  and  velocity  dependence,  and 
determining  the  relative  contributions  of  the  various  chemiluminescent 
reactions. 

As  shown  in  Section  II,  the  amount  of  radiation  emitted  per  unit 
volme  in  the  wake  by  the  three  principal  chemiluminescent  reactions, 
reactions  (1),  (2),  and  (3),  will  be  given  by 

- 17  2  3 

=  3.1  X  10  photons/cm  ‘sec 

=  1.9  X  10  +  2.6  X  10  photons /cm^- sec 

photons/cm^-sec 

where  refers  to  reaction  (i),  and  the  concentrations  are  to  be 

3 

expressed  in  particles /cm  .  The  amount  of  radiation  frcaa  reaction  (2) 
o 

emitted  above  3000  A  will  also  be  of  interest.  As  discussed  in 
Section  II,  this  is  given  by 

J2  =  0.9  X  10  ^  photons /cm^*  sec 


J  =  6.4  X  10-  Nq 
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The  total  amount  of  wake  radiation  emitted  by  these  three  reactions 
will  depend  on  the  species  concentrations  and  on  the  effective  radiating 
volume  of  the  wake.  In  order  to  bring  out  the  important  features  of 
the  problem  without  unduly  complicating  the  calculation,  we  will 
assume  that  the  axis  values  of  the  species  concentrations  hold  over 
a  cylindrical  region  surrounding  the  axis.  Bie  total  amount  of  wake 
radiation  from  reaction  (i)  will  then  be 


I. 

X 


TT  L.  J, 
1  i 


where  ir  R  is  the  effective  cross  sectional  area  of  the  radiating 

wake,  and  the  wake  length  appropriate  to  reaction  (i). 

As  shown  by  the  results  of  Lin  and  Hayes, for  reactions 

involving  N  or  0,  R  ~  3  R^.  The  effective  wake  length  will  depend  on 

the  reactants  producing  the  radiation.  For  reaction  (3),  NO  +  0  -♦  NO2  +  hu, 

the  effective  wake  length  will  be  the  length  of  the  atomic  oxygen  wake. 

This  will  be  given  by  L  =  wake  downstream  of  the  P 

point  is  moving  at  very  nearly  the  freestream  velocity.  For  reaction  (1), 

N  +  N  N2  +  hn,  and  reaction  (2),  N  +  0  -»  NO  +  hu,  the  effective  wake 

length  will  be  determined  by  the  length  of  the  atomic  nitrogen  wake. 

As  shown  in  Section  III,  this  will  depend  on  the  amount  of  turbulent 

mixing.  In  the  no  mixing  limit,  it  will  be  given  by  L  =  where 

we  use  the  values  of  Fig.  10.  In  the  limit  of  homogeneous  mixing, 

the  atomic  nitrogen  wake  length  will  be  determined  by  the  downstream 

distance  required  for  the  turbulent  core  to  expand  into  the  undissociated 

2  (1) 

flow.  This  will  be  approximately  L  ~  10 
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Using  this  model  of  the  radiating  wake,  in  Figs.  13-18  we  present 

the  results  for  a  nose  radius  of  1  ft,  using  the  two  extreme  models 

for  the  turbulent  mixing  process.  In  Fig.  13  we  show  the  radiation 

frcm  the  Lewis-Rayleigh  nitrogen  afterglow  (Ng  1st  positive  bands), 

in  Fig.  14  the  nitric  oxide  afterglow  (1^  y,  and  8  bands),  and  in 

Fig.  15  the  1K)2  continuum.  Figure  16  shows  the  radiation  frcmi  the 

o 

nitric  oxide  afterglow  above  3000  A.  In  Fig.  17  we  present  the  total 

nonequilibriim  wake  radiation  from  all  three  chemiluminescent  reactions, 

o  * 

and  in  Fig.  18  the  total  radiation  above  3000  A. 

Ihe  altitude  and  velocity  dependence  are  shown  in  the  graphs. 

The  scaling  for  body  size  will  be  fairly  simple:  In  the  no  mixing 

2 

limit  the  radiation  will  scale  as  R^;  in  the  homogeneous  mixing  limit, 

o  3 

the  NO^  radiation  will  scale  as  R  ,  the  No  and  NO  radiation  as  Rr. 

2  n  Z  ^ 

We  see  from  these  results  that  the  nature  of  the  wake  radiation 

depends  very  strongly  on  the  details  of  the  turbulent  mixing  process. 

In  the  limit  of  no  microscopic  mixing  (inviscld  random  convection)  the 

N2  1st  positive  bands,  and  the  NO  bands  are  the  dcmlnant  source 

of  radiation  at  all  altitudes,  with  the  NO2  continuum  being  several 

orders  of  magnitude  less  intense.  In  the  limit  of  Instantaneous 

hcsaogeneous  mixing,  however,  the  NO2  contlnum  is  greatly  enhanced,  and 

the  N2  and  NO  radiation  reduced,  particularly  at  higher  altitudes,  so 

that  above  about  150,000  ft  the  NO2  radiation  becoms  d<Mainant.  The  total 

amount  of  radiation  is  also  affected  by  the  turbulence;  at  the  higher 

altitudes  the  difference  between  the  two  turbulent  models  can 

^  3  3 

In  changing  units  from  photons /cm  ‘sec  to  watts /cm  ,  we  have 

assumed  average  photon  energies  of  i  N2  1st  positive  syst^,  2  ev; 

NO  Y,  and  6  systems,  5  ev;  NO2  continuvma,  2  ev. 
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Fig.  13— Wake  radiation  in  N2  first  positive  bands 
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Fig. 15— Wake  radiation  in  NO2  continuum 
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Fig.  17— Total  wake  radiation 
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Fig.  18 —  Total  wake  radiation  above  3000  A 
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be  almost  an  order  of  magnitude. 

As  was  mentioned  in  Section  III,  because  of  uncertainties  concerning 
the  detailed  nature  of  the  turbulent  process  it  is  not  possible  at 
the  present  time  to  specify  quantities  affected  by  turbulent  mixing 
with  any  greater  accuracy  than  that  represented  by  the  two  turbulent 
models  used  here.  For  that  reason,  the  differences  in  the  radiation 
predictions  of  the  two  turbulent  models,  as  shown  in  Figs.  13-18, 
represent  a  real  uncertainty  in  our  ability  to  predict  the  amount  of 
nonequillbrltua  wake  radiation.  Because  of  this  uncertainty,  it  would 
be  Inappropriate  at  present  to  attempt  a  more  detailed  estimate  of  the 
amotint  of  wake  radiation  than  was  attaapted  here.  It  would  have  been 
possible,  for  example,  to  construct  detailed  profiles  of  temperature 
and  species  concentration  in  the  wake,  obtain  detailed  profiles  of 
radiation  emission  from  these,  and  then  nximerically  integrate  over  the 
CTtire  wake  to  obtain  the  total  emission.  However,  in  the  absence 
of  new  knowledge  concerning  the  turbulent  mixing  process,  the  apparent 
high  accuracy  of  such  a  calculation  would  be  specious. 

In  Fig.  19  we  compare  the  wake  radiation  with  the  gas  cap  and 

surface  radiation.  For  the  gas  cap  radiation  we  have  used  Kivel's 

(31)  3 

tables,  together  with  an  effective  radiating  volume  of  R^/10. 

For  the  surface  radiation  we  have  assumed  an  enrailsivity  of  unity, 

2 

and  an  effective  radiating  area  of  tt  r^.  We  see  from  these  results 
that  the  gas  cap  radiation  will  be  negligible  compared  to  the  wake 
radiation  for  all  altitudes  considered  here.  ,The  surface  radiation 
will  become  important  once  the  surface  tsnperature  approaches  3000°K. 
However,  this  will  most  probably  not  happen  until  relatively  low 
altitudes,  100,000  ft  or  less,  so  that  the  wake  radiation  will  be 

The  values  shown  are  for  a  velocity  of  20,000  ft/sec.  The  results 
for  other  velocities  are  similar. 


80,000  100,000  120,000'  140,000  160,000  180,000  200,000  220,000 

Altitude  (ft) 


Fig.l9— Comparison  of  wake,  gas  cap,  and  surface  radiation 
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dominant  at  the  higher  altitudes. 

In  using  the  results  obtained  in  this  Meanorandum  it  should  be 
remembered  that  we  have  been  considering  a  pure  air  wake.  For  wakes 
containing  ablation  products,  which  are  much  more  effective  radiators 
than  air,  the  wake  radiation  would  probably  be  greatly  enhanced  over 
the  values  shown  here. 


V.  DISeUSSIQN 


In  this  Memorandum  we  have  considered  the  chemiluminescent 
reactions  that  will  be  important  sources  of  wake  radiation  under 
nonequilibrim  conditions,  we  have  discussed  the  various  aerodynamic 
and  chemical  phenomena  that  will  affect  the  radiation,  and  we  have 
made  estiiaates  of  the  amount  of  radiation  emitted.  The  most  striking 
fact  to  emerge  from  this  study  has  been  the  sensitivity  of  the 
nonequilibriinn  radiation  to  details  of  the  turbulent  mixing  process. 
The  present  study  has  been  limited  to  pure  air  wakes,  but  it  is  almost 
certain  that  the  nonequilibrium  radiation  from  contaminated  wakes 
will  be  equally  sensitive  to  turbulence.  This  dependence  of  the  wake 
radiation  upon  presently  uncertain  details  of  the  turbulent  process 
provides  another  reason,  if  any  were  needed,  for  the  importance  of 
increasing  our  knowledge  of  hypersonic  turbulence. 

Considering  the  sensitivity  of  the  nonequilibrium  wake  radiation 
to  changes  in  the  amount  of  turbulent  mixing,  it  is  possible  that 
radiation  measurements  could  be  used  to  study  the  turbulent  process. 

For  example,  if  reentry  wakes  of  nonablating  vehicles  could  be 
observed  with  sufficient  spectral  resolution  to  measure  the  relative 
intensities  of  the  N2  first  positive  system,  the  NO  bands,  and  the 
NO2  continuum,  it  is  conceivable  that  the  results  could  be  used  to 
choose  between  the  Invlscid  random  convection  model  and  the  homogeneous 
mixing  TK>del,  or  to  establish  a  new,  intermediate  model  for  the 
turbulent  process.  Observations  of  the  wake  length  of  the  N2  or  NO 
bands  could  also  be  used  for  this  purpose.  Similar  observations  of 
the  wakes  of  ablating  vehicles  could  also  be  used  to  study  the 


turbulent  process,  providing  that  the  necessary  theoretical  study  of 
nonequilibrium  radiation  from  contaminated  wakes  had  first  been 
carried  out^ 


59 


REFERENCES 


1.  Lin,  S.C.,  and  J.E.  Hayes,  A  Qnasi  One- Dimensional  Model  for 

Chemically  Reacting  Turbulent  Wakes  of  Hypersonic  Objects. 
Avco-Everett  Research  Laboratory,  Research  Report  157,  July 
1963. 

2.  Young,  R. A. ,  and  R.L.  Sharpless,  "Chemiluminescent  Reactions 

Involving  Atcanic  Oxygen  and  Nitrogen,"  J.  Chem.  Phys.  ,  Vol.  39, 
No.  4,  August  1963,  pp.  1071-1102. 

3.  Bayes,  K. D.,  and  G.B.  Kistiakowsky,  "On  the  Mechanism  of  the 

Lewis-Rayleigh  Nitrogen  Afterglow,"  J.  Chem.  Phys. ,  Vol.  32, 

No.  4,  April  1960,  pp.  992-1000. 

4.  Young,  R.A. ,  and  R.L.  Sharpless,  "Excitation  of  the  P,yj6 

Ogawa  Bands  of  Nitric  Oxide  in  the  Association  of  Atomic 
Nitrogen  and  Oxygen,"  Discussions  Faraday  Soc. .  No.  33,  1962, 
pp.  228-256. 

5.  Broida,  H.P.,  H.I.  Schiff,  and  T.M.  Sugden,  "Observations  on  the 

Chemiluminescent  Reaction  of  Nitric  Oxide  with  Atomic  Oxygen," 
Trans.  Faraday  Soc. ,  Vol.  57,  Part  2,  February  1961,  pp.  259-265. 

6.  Fontijn,  A.,  C.B.  Meyer,  and  H.I.  Schiff,  "Absolute  Quantum  Yield 

Measurements  of  the  NO-0  Reaction  and  Its  Use  as  a  Standard  for 
Chemiluminescent  Reactions,"  J.  Chem.  Phys. ,  Vol.  40,  No.  1, 
January  1964,  pp.  64-70. 

7.  Wurster,  W.H. ,  and  P.V.  Marrone,  Study  of  Infrared  Emission  in 

Heated  Air,.  Cornell  Aeronautical  Laboratory,  Report  QM-1373-A-4, 
July  1961. 

8.  Wurster,  W.H. ,  Study  of  Infrared  Ehiission  from  Hypersonic  Air 

Flows ,  Cornell  Aeronautical  Laboratory,  Report  QM-1626-A-13, 

July  1963. 

9.  Gilmore,  F.R. ,  Approximate  Radiation  Properties  of  Air  Between 

2000  and  8000 OR.  The  RAND  Corporation,  RM-3997-ARPA,  March  1964. 

10.  Polanyi,  J.C. ,  "Infrared  Chemiluminescence."  J.  Quantitative 

Spectroscopy  and  Radiative  Transfer.  Vol,  3,  No.  4,  October/ 
December  1963,  pp.  471-496. 

11.  Morgan,  J.E. ,  L.F.  Phillips,  and  H.I.  Schiff,  "Studies  of 

Vibrational ly  Excited  Nitrogen  Using  Mass  Spectrometric  and 
Calorimeter-Probe  Techniques,"  Discussions  Faraday  Soc. , 

No.  33,  1962,  pp.  118-127. 


60 


12.  Tearne,  J.D. j  S.  Georgiev,  and  R.A.  Allen,  Radiation  from  the 

Non-Equilibritim  Shock  Front,  Avco-Everett  Research  Laboratory, 
Research  Report  112,  October  1961. 

13.  Feldman,  S. ,  Hypersonic  Gas  Dynamic  Charts  for  Equilibrium  Air, 

Avco-Everett  Research  Laboratory,  Research  Report  40,  January 
1957. 

14.  Gilmore,  F.R.  ,  Equilibrium  Composition  and  Thermodynamic  Properties 

of  Air  to  24,000*^,  The  RAND  Corporation,  RM-1543,  August  1955. 

15.  Lees,  L. ,  and  T.  Kubota,  "Inviscid  Hypersonic  Flow  over  Blunt- 

Nosed  Slender  Bodies,"  J.  Aeronautical  Sciences,  Vol.  24, 

March  1957,  pp.  195-197. 

16.  Feldman,  S. ,  A  Numerical  Comparison  between  Exact  and  Approximate 

Theories  of  Hypersonic  Inyiscid  Flow  Past  Slender  Blunt-Nosed 
Bodies,  Avco-Everett  Research  Laboratory,  Research  Report  71, 

June  1959. 

17.  Goulard,  M. ,  and  R.  Goulard,  The  Aerothermodynamics  of  Reentry 

Trails,  American  Rocket  Society,  Preprint  1145-60,  Ifay  1960. 

18.  Hall,  J.G. ,  A.Q.  Eschenroeder,  P.V.  Marrone,  Inviscid  Hypersonic 

Airflows  with  Coupled  Nonequilibrium  Processes,  Cornell 
Aeronautical  Laboratory,  Report  AF-1413-A-2,  May  1962. 

19.  Geiger,  R.E. ,  "On  the  Frozen  Flow  of  a  Dissociated  Gas," 

J.  Aerospace  Sciences,  Vol.  26,  December  1959,  pp.  834-835. 

20.  Lin,  S.C.,  and  J.D.  Teare,  A  Streamtube  Approximation  for 

Calculation  of  Reaction  Rates  in  the  Inviscid  Flow  Field  of 
Hypersonic  Objects,  Avco-Everett  Research  Laboratory,  Research 
Note  223,  August  1961. 

21.  Wray,  K.L. ,  Chemical  Kinetics  of  High  Temperature  Air,  Avco- 

Everett  Research  Laboratory,  Research  Report  104,  June  1961. 

22.  Bortner,  M. H.,  Chemical  Kinetics  in  a  Reentry  Flow  Field, 

General  Electric  Missile  and  Space  Division,  Report  R63SD63, 

August  1963. 

23.  Lin,  S.C.,  and  J.D.  Teare,  Rate  of  Ionization  Behind  Shock 

Waves  in  Air,  II,  Theoretical  Interpretation,  Avco-Everett 
Research  Laboratory,  Research  Report  115,  September  1962. 

24.  Back,  R.A. ,  W.  Dutton,  C.A.  Winkler,  "The  Decay  of  Active  Nitrogen 

at  High  Temperature,"  Canadian  J.  Chem. ,  Vol.  37,  December  1959, 
pp.  2059-2063. 

25.  Harteck,  P. ,  and  R. R.  Reeves,  "Recent  Investigations  of 

Chemical  Reactions  of  Fundamental  Importance  in  the  Atmosphere," 
in  R. D.  Cadle  (ed.),  Chemical  React ions  in  the  bower  and  Upper 
Atmosphere,  Chap.  14,  Interscience  Publishers,  1961,  pp.  219-238. 


26. 


Feldraan,  S» ,  Trails  of  AxirSvngaetric  Hypersonic  Blunt  Bodies 

Flying  through  tke  Atmosphere,  Avco-Everett  Res ear eh  Laboratory, 
Research  Report  82,  Beeember  1959. 

27.  Lees,  L. ,  and  L.  Hromas,  "Turbulent  Diffusion  in  the  Wake  of  a 

Blunt-Nosed  Body  at  Hypersonic  Speeds,"  J.  Aerospace  Sciences, 
Vol.  29,  August  1962,  pp.  976-993. 

28.  Lin,  S.G. ,  Ionized  Wakes  of  Hypersonic  Objects,  Avco-Everett 

Research  Laboratory,  Research  Report  151,  June  1959. 

29.  Webb,  W.H, ,  L.  Hromas,  and  L.  Lees,  "Hypersonic  Wake  Transition," 

AIM  J.  %  Vol.  1,  No.  3,  March  1963,  pp.  719-721. 

30.  Zeiberg,  S.L. ,  and  G.D.  Bleich,  "Finite  Difference  Calculation 

of  Hypersonic  Wakes,"  paper  presented  at  AIM  Conference  on 
Physics  of  Entry  into  Planetary  Atmospheres,  Cambridge, 
Massachusetts,  August  26-28,  1963. 

31.  Kivel,  B. ,  Radiation  from  Hot  Air  and  Stagnation  Heating, 

Avco-Everett  Research  Laboratory,  Research  Report  79,  October 
1959. 


